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It is widely believed that when a molecule with thiol (S-H) end groups bridges a pair of gold 
electrodes, the S atoms bond to the gold and the thiol H atoms detach from the molecule. However, 
little is known regarding the details of this process, its time scale, and whether molecules with and 
without thiol hydrogen atoms can coexist in molecular junctions. Here we explore theoretically how 
inelastic tunneling spectroscopy (IETS) can shed light on these issues. We present calculations of 
the geometries, low bias conductances and IETS of propanedithiol and propanedithiolate molecular 
junctions with gold electrodes. We show that IETS can distinguish between junctions with molecules 
having no, one or two thiol hydrogen atoms. We find that in most cases the single- molecule junctions 
in the IETS experiment of Hihath et al. [Nano Lett. 8, 1673 (2008)] had no thiol H atoms, but that a 
molecule with a single thiol H atom may have bridged their junction occasionally. We also consider 
the evolution of the IETS spectrum as a gold STM tip approaches the intact S-H group at the end of 
a molecule bound at its other end to a second electrode. We predict the frequency of a vibrational 
mode of the thiol H atom to increase by a factor ^ 2 as the gap between the tip and molecule 
narrows. Therefore, IETS should be able to track the approach of the tip towards the thiol group 
of the molecule and detect the detachment of the thiol H atom from the molecule when it occurs. 

PACS numbers: 81.07.Nb, 72.10.Di, 73.63.Rt, 85.65.+h 



I. INTRODUCTION 

Electrical conduction in molecular junctions consist- 
ing of a single organic molecule chemically bonded to a 
pair of metal contacts has been studied experimentally 
and theoretically for more than a decade. [1] However, di- 
rect experimental observation of the atomic scale struc- 
tures of the molecule-electrode interfaces in such devices 
is not feasible at the present time and the determina- 
tion of the nature of these interfaces continues to be an 
important goal in this field. [1] The most studied single 
molecular junctions are those in which the molecule is 
terminated at its ends by thiol (S-H) groups and bonds 
to a pair of gold electrodes via the sulfur atoms. [1] It 
is widely believed that when such a molecule bonds to 
gold electrodes, the hydrogen atoms of the thiol groups 
are displaced by the gold and leave the vicinity of the 
molecule so that the thiol groups become thiolates.[l, 2] 
However, ab initio calculations [6, 8, 9] have indicated 
that the energies of such molecular junctions with and 
without hydrogen atoms attached to the sulfur atoms are 
very similar, and that junctions with intact S-H groups 
might even have the lower energies [8, 9], suggesting that 
some degree of coexistence of species with intact and dis- 
sociated S-H groups bound to gold electrodes may be 
possible. [6-9] Furthermore whether the cleavage of the 
S-H bond occurs on the time scales on which single- 
molecular junctions are formed (and broken) in statistical 
STM break junction experiments [3-5] (an experimental 
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technique widely used today in studies of single-molecule 
nanoelectronics [1]) has still not been established. It is 
therefore of interest to develop a methodology that can 
determine unambiguously whether a thiol hydrogen atom 
is present in the gold-molecule- gold junction at a given 
time and can provide some information regarding the 
atomic structure of the thiol-gold region during the pro- 
cess by which the chemical bond between the molecule 
and the gold electrode forms. In this paper we explore 
this issue theoretically and propose inelastic tunneling 
spectroscopy (IETS) as an attractive candidate for such 
a methodology. 

In single-molecule IETS experiments a bias voltage is 
applied across a molecular junction at low temperatures 
and a conductance step is observed in the current- volt age 
characteristic of the system when the applied voltage 
increases beyond the threshold value for the excitation 
of a vibrational quantum (commonly referred to as a 
"phonon") in the junction. Many experimental IETS 
studies of molecular junctions have been carried out [10- 
40] and theoretical work [41-77] has accounted for vari- 
ous aspects of the IETS data. These IETS experiments 
together with the theoretical work demonstrated conclu- 
sively that particular molecular species are involved in 
electrical conduction through metal-molecule- metal junc- 
tions. It has also been shown that IETS can be sensitive 
to the molecular conformation [29, 38, 74, 78], the ori- 
entation of molecules relative to the electrodes [56, 77] 
and the gold-sulfur atomic bonding geometries at the 
molecule- metal interfaces [75, 76, 78]. 

However, none of the above studies examined the effect 
that thiol hydrogen atoms, if present, might have on the 
inelastic tunneling spectra of molecular junctions. We do 
this in the present article, considering as a specific exam- 
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pie molecular junctions with gold electrodes bridged by 
molecules derived from 1,3-propanedithiol (PDT) that re- 
tain one or both thiol hydrogen atoms. We also compare 
our results for these systems with our previous findings 
[75, 78] for 1,3-propanedithiolate molecules (with no thiol 
hydrogen atoms) bridging gold electrodes and with the 
experimental statistical STM break junction IETS data 
of Hihath et al [26]. We find IETS to be able to dis- 
tinguish molecules bridging gold electrodes with no hy- 
drogen atoms bound to the sulfur atoms of the molecule 
from molecules with a hydrogen atom bound to one of the 
sulfur atoms, and from molecules with hydrogen atoms 
bound to both sulfur atoms. Our results indicate that 
junctions with no hydrogen atom bound to either sulfur 
atom of the bridging molecule were the systems probed 
in most cases by the IETS measurements of Hihath et 
al. [26]. We also find the data of Hihath et al. to show no 
evidence of junctions with hydrogen atoms bound to both 
sulfur atoms of the bridging molecule but to be consistent 
with transport through molecules with a hydrogen atom 
bound to only one of the sulfur atoms being observed 
occasionally. We also find that a vibrational mode asso- 
ciated with the thiol hydrogen atom has a strong IETS 
signature and that its frequency is sensitive to the dis- 
tance between the thiol sulfur atom and the gold elec- 
trode to which that sulfur atom bonds. We therefore 
propose that by tracking the frequency of this vibrational 
mode in a statistical STM break junction experiment it 
should be possible to monitor the closing of the gap be- 
tween the thiol-terminated end of the molecule and a gold 
electrode in real time and to detect the detachment of the 
thiol hydrogen atom from the molecule. 



The remainder of this paper is organized as follows: 
In Section II, we outline how we calculate the relaxed 
molecular junction geometries, vibrational modes, low 
bias elastic tunneling conductances, and the conductance 
step heights associated with inelastic electron scattering 
in the molecular junctions. In Section III we present our 
results. In Section III A We discuss the conformations of 
PDT molecular junctions with different numbers of sulfur 
atoms having attached hydrogen atoms that we obtained 
by density functional theory-based unconstrained relax- 
ations of the containing one or two molecules. Our cal- 
culated low bias elastic conductances of these junctions 
are discussed in Section III B. In Section III C, we present 
the results of our calculations of the inelastic tunneling 
spectra of these junctions and compare them with the ex- 
perimental measurements of Hihath et al. [26] In Section 
III D we present our results regarding the evolution of the 
inelastic tunneling spectrum of a PDT molecule one end 
of which is bound to a gold electrode while an intact S-H 
group at the other end is approached by a gold STM tip. 
We present our conclusions and discuss the experimental 
implications of our results in Section IV. 



II. THEORY 

The theoretical approach that we adopt here is sim- 
ilar to that used in Refs. 75 and 78. In the work re- 
ported in those papers [75, 78] this approach was very 
successful in accounting for the relevant experimental 
IETS data [26] and identifying for the first time the ex- 
perimentally realized gold-sulfur bonding geometries in 
gold-propanedithiolate-gold molecular junctions. An in 
depth discussion of the underlying theory and the rel- 
evant mathematical derivations have been presented in 
Ref. 78. We therefore present only a brief summary of 
the theoretical approach that we use here and refer the 
reader to Ref. 78 for further details. 

In this work we focus on inelastic tunneling processes 
that involve vibrational modes with large amplitudes of 
vibration of the hydrogen and sulfur atoms that are lo- 
cated at the thiol-gold interface. Therefore, it is nec- 
essary to calculate accurate relaxed geometries of both 
the molecule and nearby atoms of the gold electrodes as 
well as accurate atomic displacements from equilibrium 
for these atoms in the vibrational modes. We do this by 
modeling the molecular junctions as extended molecules 
that consist of the molecule itself and two clusters of gold 
atoms representing the pair of electrodes to which the 
molecule binds. Unlike many previous theoretical stud- 
ies that calculated the relaxed geometries and vibrational 
modes of only the molecule itself, keeping the positions of 
the electrode atoms frozen in the geometry of a bulk gold 
crystal, we calculate the relaxed geometries and vibra- 
tional modes of the entire extended molecules including 
both the molecule and and the gold atoms. The relaxed 
geometries and vibrational normal modes of these ex- 
tended molecules were computed using density functional 
theory[79, 80] as implemented in the GAUSSIAN 09 [81] 
package with the B3PW91 density functional [82] and the 
Lanl2DZ pseudo-potentials and basis set. [81] Our calcu- 
lations were repeated for gold clusters of different sizes 
with up to 14 gold atoms per cluster. These gold cluster 
sizes were found to be sufficient to achieve convergence of 
our results for the vibrational mode frequencies and in- 
elastic tunneling intensities with increasing cluster size. 
A physical reason for this rapid convergence with increas- 
ing cluster size is that the frequencies of the vibrational 
modes that we study are higher than the frequency of 
the highest phonon mode of bulk gold and for this rea- 
son the vibrational modes that we study have a strongly 
evanescent character in the gold clusters. That is, their 
vibrational amplitudes in the gold clusters decay very 
rapidly with increasing distance from the site to which 
the molecule bonds and are negligible on the furthest 
gold atoms from the molecule for the larger gold clusters 
that we consider. In this article we consider fully relaxed 
extended molecules with no geometrical constraints im- 
posed, and also extended molecules subject to an exter- 
nally imposed strain that we model by keeping the dis- 
tance between two outer atoms of the gold clusters fixed 
during the relaxation of the extended molecule. We note 
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that the experimental inelastic tunneling spectra of Hi- 
hath et al. [26] were taken by sweeping the bias applied 
across the molecular junction while keeping the STM tip 
stationary relative to the substrate. 

Our calculations of the low bias elastic conductances of 
the molecular junctions combine the results of the above 
ab initio calculations of the extended molecule geometries 
with Landauer transport theory, [1] Green's function tech- 
niques, the solution of the Lippmann-Schwinger equation 
and the semi-empirical extended Hiickel model with the 
parameters of Ammeter et al. [83] to estimate the elec- 
tronic structures for these geometries. As is discussed 
in detail in Refs.l and 85, this methodology involves 
no fitting to any experimental data relating to trans- 
port in molecular junctions. It is known to yield low 
bias conductances in good agreement with experiments 
for propanedithiolate bridging gold electrodes [75, 78] as 
well as other molecules thiol-bonded to gold electrodes 
[1, 84, 85]. We evaluate the zero bias tunneling conduc- 
tances from the Landauer formula [1] 

. . 1 1 u i 

with the quantum unit of conductance go = 2e 2 /h, after 
evaluating the elastic transmission amplitudes tj\ from 
the solution of the Lippmann-Schwinger equation. In 
the transmission amplitude tji, i is the electronic state of 
a carrier that is coming from the left lead, and j is the 
electronic state of a carrier that has been transmitted to 
the right lead, vj and Vi are the corresponding electron 
velocities. The coupling of the extended molecule to the 
electron reservoirs was treated as in previous work [75, 
78, 85-94] by attaching a large number of semi-infinite 
quasi-one-dimensional ideal leads to the valence orbitals 
of the outer gold atoms of the extended molecule. 

We calculate the IETS intensities (i.e., the conduc- 
tance step heights) Sg a associated with the emission of 
phonons of mode a perturbatively, in the spirit of an ap- 
proach proposed by Troisi et al. [46] who transformed the 
problem of calculating IETS intensities into an elastic 
scattering problem. However, unlike Troisi et al. [46], we 
formulate IETS intensities explicitly in terms of elastic 
electron transmission amplitudes tj\ through the molec- 
ular junction. As is shown in Ref. 78, we find 

e 2 ^^| ^({Ad„ Q })-^({0» 2 
og a = - lim > — — , (2) 

at low temperatures. Eq. 2 states that the scattering am- 
plitude for inelastic transmission of an electron through 
the molecular junction is proportional to the change in 
the elastic amplitude for transmission through the junc- 
tion when its atoms are displaced from their equilibrium 
positions due to the excitation of vibrational mode a. 
Here ^-({0}) is the elastic electron transmission ampli- 
tude through the molecular junction in its equilibrium 
geometry from state i with velocity V{ in the electron 



source to state j with velocity vj in the electron drain. 
d na represents the displacements from their equilibrium 
positions of the atoms n of the extended molecule in nor- 
mal mode a normalized so that m n d* a / -d na = 5 a ^ a 
where m n is the mass of atom n. u a is the frequency 
of mode a. tj\({Ad na }) is the elastic electron transmis- 
sion amplitude through the molecular junction with each 
atom n displaced from its equilibrium position by Ad nai 
where A is a small parameter. 

We evaluate tj[ in Eq. (2) numerically to find the 
heights Sg a of the conductance steps that arise from in- 
elastic tunneling processes due to emission of phonons of 
vibrational mode a. The calculations are carried out at 
the zero bias Fermi energy since the values of the bias 
at which the inelastic transmission occurs in the experi- 
ments of Hihath et al. [26] are low. 

The validity of our methodology rests largely on the 
fact that we rely on ab initio density functional theory 
to calculate the frequencies and atomic displacements 
from equilibrium for the vibrational normal modes that 
give rise to the features of interest in the inelastic tun- 
neling spectra. Such calculations are accurate because 
in the Born-Oppenheimer approximation they consti- 
tute electronic ground state total energy calculations for 
which density functional theory is appropriate [1, 79, 80] 
and for which the density functional model that we 
use [81, 82] has been optimized. Our key predictions con- 
cern these accurately calculated frequencies (or the asso- 
ciated phonon energies) for the modes with the largest 
intensities in the inelastic tunneling spectra, and the de- 
pendence of these frequencies on the details of the confor- 
mation of the molecular junction. We also calculate the 
junction conformations using density functional theory. 
We rely on the extended Hiickel model and perturba- 
tion theory only to determine which of the vibrational 
modes have the largest inelastic tunneling intensities. It 
should also be noted that perturbative calculations of 
this type are appropriate provided that the transmitted 
electrons have energies far from any elastic resonance 
and if the amplitudes of the excited vibrational modes 
are small. These conditions are expected to be satisfied 
at low bias voltages and cryogenic temperatures such as 
those in the experiment of Hihath et al. [26] Also, the ex- 
tended Hiickel model that we use has been found to yield 
results in reasonably good agreement with experiment 
for elastic[l, 75, 78, 84, 85, 95, 96] and inelastic[75, 78] 
transport in a variety of molecular wires thiol-bonded to 
gold electrodes. 



III. RESULTS 

A. Conformations of extended molecules obtained 
by unconstrained relaxations 

We considered extended molecules where the molecule 
bonds via its sulfur atoms to finite clusters of gold atoms, 
and relaxed these entire structures using density func- 



Relaxed Structures 




FIG. 1. (Color online.) Examples of relaxed structures 
of trans-PDT molecules bonded to gold clusters in the top 
geometry. [Ill] (a): A PDT molecule with no hydrogen atom 
attached to either sulfur atom, (b): A PDT molecule with 
a hydrogen atom attached to one of the sulfur atoms, (c): 
A PDT molecule with a hydrogen atom attached to each of 
the two sulfur atoms, (d): Two PDT molecules bridging the 
gold clusters with no hydrogen atom attached to any sulfur 
atom, (e): Two PDT molecules bridging the gold clusters 
with a hydrogen atom attached to one of the sulfur atoms of 
each molecule. The sulfur atom of the molecule on the right 
has not chemisorbed to the lower gold cluster while each of 
the other three sulfur atoms in the system has chemisorbed 
to a gold cluster, (f): Two PDT molecules bridging the gold 
clusters with a hydrogen atom attached to each of the two 
sulfur atoms of each molecule. Carbon, hydrogen, sulfur and 
gold atoms are black, blue, yellow and amber, respectively. 



tional theory-based calculations. [79-82] We carried out 
the density functional geometry relaxations starting from 
initial trial geometries of the PDT molecule with a hy- 
drogen atom attached to one, both or neither of the two 
sulfur atoms of the molecule. 

We were able to obtain relaxed structures in which 
the hydrogen atom(s) remaining attached to one or both 
of the sulfur atom(s) for bonding geometries in which 
the sulfur with the attached hydrogen atom bonds to 
one gold atom (referred to as "top site" bonding geome- 



TABLE I. (Color online.) Calculated low bias elastic conduc- 
tances, and mode I phonon energies and IETS intensities for 
the relaxed structures of one and two molecules with different 
numbers of thiol hydrogen atoms bridging gold atomic clus- 
ters. The numerical values labelled (a) -(f) are for the molec- 
ular junctions shown in Fig. 1 (a) - (f), respectively. The 
experimentally measured low bias conductance values of Hi- 
hath et al. for single-molecule PDT junctions were ~ 0.006go 
±0.002. [26] go 2e 2 /h. 



Configuration 


Low Bias 
Conductance 

(so) 


Phonon 
Energy 
of Mode-I 

(meV) 


IETS 
Intensity 
of Mode-I 

(10" 6 5o) 


a) 


0.0013 


44.88 


10.03 


b) 


0.0017 


40.17 


4.04 


c) 


0.0021 


36.7 


6.20 


d) 


0.0029 


45.13 


10.14 


e) 


0.0019 


40.67 


3.02 


f) 


0.0041 


36.95 


9.04 



tries) but not for geometries in which the sulfur bonds 
to two or three gold atoms. Therefore in the remainder 
of this paper we will consider only top site bonded ge- 
ometries; the two other types of bonding geometries for 
PDT molecules with no attached thiol hydrogen atoms 
have been discussed in Refs 75 and 78. 

Representative examples of relaxed PDT extended 
molecules with the trans conformation are shown in 
Fig. 1. In Fig. 1, all S atoms bond to the gold clusters 
in the top site geometry. In all of the examples shown 
in Fig. 1 each gold cluster contains 14 atoms and the re- 
laxations were performed with no constraints applied to 
the extended molecule geometries. In Fig. 1 (a), (b) and 
(c) a single PDT molecule with zero, one and two thiol 
hydrogen atoms bridges the gold clusters. In Fig. 1 (d), 
(e) and (f) two PDT molecules, each with zero, one and 
two thiol hydrogen atoms, bridge the gold clusters. 

Our findings regarding the energetics of the structures 
shown in Fig. 1 are qualitatively similar to those of the 
previous studies [6-9] that have indicated that the ener- 
gies of molecular junctions with and without hydrogen 
atoms attached to the sulfur atoms are very similar, and 
that junctions with intact S-H groups might even have 
the lower energies. For example, our density functional 
theory-based calculations yielded an energy for the struc- 
ture in Fig. 1 (a) (that has no thiol H atoms) that, to- 
gether with the energy of a free H2 molecule, is higher 
than the energy of the structure in Fig. 1 (c) (that has two 
thiol H atoms) by 0.27eV. Similarly, we found an energy 
for the structure in Fig. 1 (d) (that has no thiol H atoms) 
that, together with the energy of two free H2 molecules, 
is higher than that of the structure in Fig. 1 (f) (that has 
four thiol H atoms) by 0.04eV. This suggests that a sys- 
tem with a dithiol molecule bridging gold electrodes may 
be slightly more stable energetically than that formed if 
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the molecule dissociates into the dithiolate and a free H 2 
molecule. However, this ignores entropy considerations, 
which favor the latter system and furthermore it is not 
known whether the H atoms that detach from the thiol 
groups form H 2 molecules or other, possibly more stable, 
moieties. It is also worth noting that the geometries of 
the gold clusters in Fig. 1 (a) and (c) are not identical, 
as is also the case for Fig. 1 (d) and (f), and that con- 
clusions drawn from the above comparisons between the 
energies of these structures should be approached with 
some caution for this reason as well. 



B. Elastic conductances in the limit of low bias 

Our calculated zero bias conductance values for the 
molecular junctions shown in Fig. 1 are presented in in 
Table I where rows (a) - (f) correspond to structures (a) 
- (f) of Fig. 1, respectively. 

As is seen in Table I, the calculated zero bias conduc- 
tance values are g=0. 0013go, 0.0017go and 0.002 lgo for 
top-site bonded single- molecule junctions with no thiol 
hydrogen atoms, one thiol hydrogen atom, and two thiol 
hydrogen atoms, respectively. All of these theoretical val- 
ues agree with the experimental value 0.006 ± 0.002#o of 
Hihath et al. [26] to a degree that is typical of the exper- 
imental and theoretical literature [97] for molecules thiol 
bonded to gold electrodes. These calculated elastic con- 
ductance values show an increasing trend with increasing 
number of thiol hydrogen atoms in the molecular junc- 
tion. However, since the calculated values differ by very 
modest amounts (less than a factor of two), determining 
experimentally whether or not hydrogen atoms remain 
bound to the sulfur atom(s) in experimental realizations 
of gold-PDT-gold molecular junctions on the basis of low 
bias conductance measurements alone is expected to be 
difficult or impossible. As shown in Section IIIC, com- 
paring our predictions for the inelastic tunneling spectra 
with experiment [26] provides strong evidence that PDT 
molecular junctions with no thiol hydrogen atoms pre- 
dominated in the experiments of Hihath et al. [26]. 

In the majority of cases the calculated zero bias con- 
ductance values in Table I for molecular junctions with 
pairs of molecules bridging the gold contacts in parallel 
are roughly a factor of 2 larger than the conductances 
of the corresponding single molecule junctions, as might 
be expected for pairs of molecules [78, 98-106] or other 
quantum conductors [107] connecting electron reservoirs. 
An exception is the structure in Fig. 1(e) for which the 
calculated low bias conductance is close to that of the 
corresponding single-molecule junction in Fig. 1(b). We 
attribute this different behavior to the molecule on the 
right in Fig. 1(e) not having formed a chemical bond be- 
tween the sulfur atom with the attached hydrogen atom 
and the nearby gold cluster, resulting in poor conduction 
through that molecule. 



C. Inelastic Tunneling Spectroscopy of Thiol 
Groups in Unconstrained Molecular Junctions 

We will begin by considering the IETS of extended 
molecules that have been relaxed without any geometri- 
cal constraints being imposed during the relaxation pro- 
cess as discussed in Section III A. Since our objective is 
to understand how the presence of thiol hydrogen atoms 
should affect the inelastic tunneling spectra of molecular 
junctions we will focus on the IETS signatures of those 
vibrational modes that have the strongest amplitudes of 
vibration on S atoms or thiol H atoms and the strongest 
IETS intensities associated with those modes. We find 
the relevant modes to fall within the energy range from 
25 to 60 meV. 

Representative examples of these vibrational modes 
and their energies are shown in Fig. 2(a) for PDT sulfur 
atoms bound to gold clusters in top site conformations (as 
discussed in Section III A) with a hydrogen atom bound 
to one or both sulfur atoms. The modes are classified ac- 
cording to the nature their most prominent atomic vibra- 
tions. In modes I, II and III the largest vibrational am- 
plitudes are those of the sulfur atoms. In mode I and II 
the two sulfur atoms move in antiphase and in phase, re- 
spectively, approximately along the axis of the molecule. 
(These modes are sometimes referred to as the symmet- 
ric and antisymmetric stretch modes, respectively.) In 
mode III the sulfur atoms move in phase, in directions 
not aligned with the molecular axis. In mode IV the most 
prominent atomic vibrations are those of the thiol hydro- 
gen atoms whose motion is roughly torsional relative to 
the direction of the nearest sulfur-carbon bond. 

The calculated IETS intensities and energies for these 
modes are shown in Fig. 2(b) together with the experi- 
mental IETS histogram of Hihath et al. [26]. Our cor- 
responding results for top site bonded PDT wires with 
no hydrogen atoms attached to the sulfur [75, 78] are 
also shown in Fig. 2(b) for comparison. Notice that 
the calculated IETS intensities corresponding to mode 
II in Fig. 2(b) are much smaller than those for the other 
modes, due to an approximate symmetry that results in 
a partial cancellation between the terms tj\ in Eq. (2) for 
mode II, as has been discussed in Refs 55, 75, and 78 for 
the case when no thiol hydrogen atoms are present. 

As can be seen in Fig. 2(b), the effect of attaching H 
atoms to the molecular S atoms on the frequency of vi- 
brational mode I (the high intensity mode that has been 
shown [75, 78] to play the dominant role in the experi- 
ment of Hihath et al. [26]) is quite dramatic: For the top 
site bound PDT junction with no H atom bound to either 
S atom, this mode occurs at 45.5 meV in Fig. 2(b) and 
is responsible for the strongest peak visible in the experi- 
mental histogram [26]. However, for the three examples of 
molecular wires each with one hydrogen atom attached to 
one sulfur atom only, mode I can be seen around 40 meV 
in Fig. 2(b), where some counts were recorded in the 
experiment [26]. The calculated IETS intensity for mode 
I for these structures is somewhat weaker than that for 
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IV 



One S is with H 



40.2 meV 




48.7 meV 






Both S are with H 




36.6 meV 




45.5 meV 





55.9 meV 




(a) 



Wavenumber(cm" 1 ) 
220 240 260 280 300 320 340 360 380 400 420 440 460 480 : . 




26 28 30 



34 36 38 40 42 44 
Energy(meV) 



54 56 58 60 



(b) 



FIG. 2. (Color online.) (a) Calculated vibrational modes 
in the energy range from 25 to 60 meV for PDT bridg- 
ing gold nano-clusters with hydrogen atoms bound to one 
or both sulfur atom(s), the sulfur atoms being bonded to 
gold in top-site geometries. These structures were relaxed 
with no constraints imposed on their geometries. Carbon, 
hydrogen, sulfur and gold atoms are black, blue, yellow and 
amber, respectively. [Ill] Arrows show un-normalized atomic 
displacements; the heavier arrows indicate the motion of the 
sulfur atoms, (b) Calculated IETS intensities (colored) vs. 
calculated vibrational mode energies for PDT molecules link- 
ing pairs of gold clusters with between 12 and 14 Au atoms in 
each cluster. Results are shown for both sulfur atoms bond- 
ing to the gold in top-site geometries and having no hydro- 
gen atom attached to either sulfur atom, having a hydrogen 
atom attached to one sulfur atom only, and having a hydrogen 
atom attached to each of the sulfur atoms. The experimental 
IETS phonon mode histogram of Hihath et al. [26] is shown 
in (darker, lighter) grey for (positive, negative) bias voltages. 
The energy ranges in which modes of types I, II, III, and IV 
in (a) occur are indicated by arrows. 



the junctions with no thiol hydrogen atoms and is very 
similar to that of mode III near 33 meV where very few 
counts were recorded in the experiment. This suggests 
while a small fraction of the PDT molecular junctions 
in the experiment [26] may have had a hydrogen atom at- 
tached to one of the sulfur atoms, most of the counts near 
40 meV in the experimental histogram [26] were probably 
not due to molecules with a hydrogen atom attached to 
one of the sulfur atoms. That is, PDT molecules with 
no thiol hydrogen atoms and with the sulfur atoms each 
bonding to two atoms of the gold electrodes (the pure 
bridge geometry) are the likely explanation of most of the 
experimental counts near 40 meV since these structures 
have strong calculated IETS intensities in that energy 
range as has been discussed in Ref.78. 

Notice that no counts are seen in the experimen- 
tal histogram [26] (Fig. 2(b)) for phonon energies near 
36.7 meV in Fig. 2(b), where we calculate mode I for the 
top site bound PDT junction with a H atom bound to 
each of the S atoms. This indicates that molecules with 
two thiol hydrogen atoms did not contribute significantly 
to the observed inelastic tunneling spectra of Hihath et 
al. [26]. 

The top site bonded molecular wires with a H atom 
bound to one or both of their S atoms all exhibit a vi- 
brational mode (mode IV in Fig. 2(a)) that involves pri- 
marily the motion of those hydrogen atoms. We find this 
mode to be located in the range 53 - 56 meV. The exper- 
imental histogram [26] shows a few counts in this range 
and its vicinity. Thus the weak feature in the experi- 
mental histogram in this range might possibly be due to 
a small minority of the molecular junctions in the exper- 
iment of Hihath et al. [26] having a hydrogen atom at- 
tached to a sulfur atom, although gauche conformations 
of the PDT molecule with no thiol hydrogen atoms may 
also have contributed to the observed inelastic tunnel- 
ing spectra in that energy range, as has been discussed 
in Ref.78. Further IETS experiments with greater sen- 
sitivity to weak features in the inelastic spectra will be 
required to discriminate definitively between these possi- 
bilities. 

In the preceding discussion we have focussed on cases 
where a single molecule connects the two electrodes. 
However in statistical STM break junction experiments 
such as that of Hihath et al. [26] the electrodes can also 
be bridged simultaneously by two (or more) molecules in 
parallel. Here we report theoretical results for the in- 
elastic contributions to the conductances of some exam- 
ples of systems of this kind. We show a pair of top site 
bonded molecules having no hydrogen atom attached to 
either sulfur atom in Fig 1(d), a pair of top site bonded 
molecules having a hydrogen atom attached to one sulfur 
atom of each molecule in Fig 1(e), and in Fig 1(f) a pair 
of top site bonded molecules having a hydrogen atom 
attached to each sulfur atom of each molecule. The in- 
termolecular distances denned as the shortest distances 
between hydrogen atoms belonging to the backbones of 
the two molecules are 4.48, 3.21, and 2.28 Angstroms re- 
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FIG. 3. (Color online.) Calculated inelastic tunneling spec- 
tra for one and two PDT molecules with different numbers of 
thiol hydrogen atoms bridging pairs of 14 atom gold clusters 
in the fully relaxed geometries depicted in Fig. 1. Labels (a) 
- (f) denote the structures (a) - (f) of Fig. 1. The atomic 
motions in modes I - IV are similar to those shown for modes 
I - IV of the single molecules in Fig. 2(a). 



spectively. The calculated inelastic tunneling spectra for 
these structures are shown together with those of the cor- 
responding single molecule junctions ( Fig l(a),(b) and 
(c) respectively) in Fig. 3. As can be seen in Fig. 3, the 
calculated inelastic tunneling spectra for these pairs of 
molecular wires are very similar to the spectra of the in- 
dividual wires apart from small frequency shifts and split- 
tings due to the weak interactions between the molecules. 
It is worth noting that although one of the molecules in 
Fig 1(e) is poorly bonded to one of the gold electrodes 
and contributes little to the low bias elastic conductance 
of the system (as has been discussed in Section IIIB) 
there is still a significant splitting of the mode I fre- 
quency for this two-molecule junction in Fig. 3 due to 
the coupling between the vibrational modes of the two 
molecules. 



D. Inelastic Tunneling Spectroscopy of Bond 
Formation Between a Gold Electrode and a 
Thiol- Terminated Molecule 

The preceding results have been for extended molecules 
that were relaxed with no geometrical constraints. Thus 
they can be regarded as pertaining to molecular junc- 
tions that have reached a local energy minimum, having 
formed chemical bonds between the molecule (s) and the 
gold electrodes. In this Section we will consider a differ- 
ent situation, namely, molecular junctions in the process 
of bond formation. More specifically, we will examine the 
evolution of the inelastic tunneling spectrum as a PDT 
molecule that has bonded to one gold electrode via a sul- 
fur atom is approached by second gold electrode (that 
may be an STM tip) so that chemical bonding between 
the molecule and the second electrode is in the process 



of taking place. 

We assume that the molecule has bonded at one of its 
ends to the first gold electrode via a sulfur atom with no 
attached hydrogen atom and that the second electrode 
approaches the thiol (S-H) group at the other end of the 
molecule. We modeled this situation by calculating the 
structures of extended molecules that are relaxed subject 
to a geometrical constraint. The role of the electrodes 
is again played by 14 atom gold clusters. We carried 
out the density functional theory-based relaxations start- 
ing from initial geometries in which the first gold cluster 
was located a typical bonding distance (about 2. 4 A) from 
the sulfur atom of the molecule with no attached hydro- 
gen atom, while the second gold cluster was located at a 
larger distance from the molecule's S-H group. During 
the relaxation process the distance I between two gold 
atoms (the furthest gold atom from the molecule in each 
of the two gold clusters) was kept fixed. By choosing dif- 
ferent values of I we obtained relaxed extended molecule 
geometries with differing values of the distance between 
the molecule's S-H group and the nearest gold cluster. 
Some representative examples of relaxed geometries ob- 
tained in this way as the gold cluster representing the 
STM tip approaches the molecule's thiol group are shown 
in Fig. 4(a) I - IV, where the gold-sulfur distances in each 
junction are also shown. A similar junction but with no 
thiol hydrogen atoms and relaxed without geometrical 
constraints is shown in Fig. 4(a) V for comparison. No- 
tice how the orientation of the PDT molecule changes 
from Fig. 4(a) I to Fig. 4(a) IV so as to assume an en- 
ergetically favorable bonding configuration in Fig. 4(a) 
where the gold-sulfur bond has fully formed. 

The calculated inelastic tunneling spectra for a series 
of structures obtained in this way are shown in Fig. 4(b). 
The S-Au distance for the sulfur atom with the attached 
hydrogen atom for each structure is shown in the legend 
together with the calculated low bias elastic conductance 
of the junction and the color used to represent the in- 
elastic tunneling spectrum for that structure. The ex- 
perimental IETS histogram of Hihath et al. [26] is also 
included in Fig. 4(b) but in this figure it serves as a ref- 
erence only since, as we have discussed in Section III C, it 
is plausible that most of the counts in the histogram were 
obtained from junctions with no hydrogen atom attached 
to either sulfur atom of the molecule being probed. The 
shaded pink areas indicate the energy ranges in which we 
find modes I, II and III (as defined in Fig. 2(a)) to occur 
for these molecular junctions. 

The most striking feature of Fig. 4(b) is the behav- 
ior of mode IV that (as in Fig. 2(a)) is characterized 
by a strong vibrational amplitude of the thiol hydrogen 
atom, and is marked by the chartreuse colored arrows in 
Fig. 4(b). The energy of this mode is ~ 26 meV (marked 
by the black arrow) when there is no gold cluster in the 
vicinity of the S-H group, i.e., the STM tip is at infinity. 
As the gold cluster approaches the S-H group the mode 
IV energy increases dramatically, reaching a value of ~ 55 
meV when the gold-sulfur distance reaches the value of 
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FIG. 4. (Color online.) (a) I to IV: Calculated geometries of PDT molecular junctions with a H atom bound to one S atom 
of the molecule. C, H, S and Au atoms are black, grey, yellow and amber, respectively. In the relaxations the separation of 
the atoms of the two Au clusters furthest from the molecule was held fixed. Au-S distances are shown. V: Molecular junction 
relaxed without constraints and with no H atom attached to either S atom, (b) Theoretical inelastic tunneling spectra of the 
PDT wires, some of which are shown in (a), color coded according to the Au-S distance for the S atom with an attached H 
atom. The calculated low bias elastic conductance of each structure is shown in the legend. Chartreuse arrows indicate mode 
IV. Black arrow indicates the energy of mode IV if there is no Au cluster near the S-H group. Pink areas indicate energy 
ranges in which modes I, II and III occur. The experimental phonon mode histogram (grey) of Hihath et al. [26] is shown as 
a reference, (c) Au-S and thiol H vibrational modes are plotted vs. Au-S distance; Au-H distances are also shown. Mode IV 
points are circled in black. The solid red curve is a fit of the mode IV frequencies to an inverse cube dependence on the Au-S 
distance; see text. 



2. 65 A at which the chemical bond between the gold and 
sulfur atoms in the presence of the thiol hydrogen atom 
has fully formed. This doubling of the mode IV energy 
as the gold STM tip approaches the molecule is a large 
effect and therefore by tracking the energy of mode IV it 
should be possible to monitor the approach of an STM 
tip to a molecule adsorbed on a surface experimentally. 
The disappearance of mode IV from the inelastic tunnel- 
ing spectrum and the simultaneous jump in the energy 
of mode I from ~ 40 meV to ~ 46 would then signal the 
ultimate detachment of the thiol hydrogen atom from the 
sulfur atom in such an experiment. 



In Fig. 4(c) we plot the wavenumbers of modes I - IV 
vs the Au - S distance as the distance between the gold 
cluster and the molecule is varied. The mode IV results 
are circled in black for clarity. While the energies of all of 
the vibrational modes plotted show some variation with 
the Au - S distance, that of mode IV is by far the largest. 
As can be seen in Fig. 4(c), this variation of mode IV fits 
well to the empirical form n = 208 + 486(x — 1.35) ~ 3 , the 
red curve in Fig. 4(c). Here n is the mode IV wavenum- 
ber in cm -1 and x is the Au - S distance in Angstoms. 
A similar (bond length to the negative third power) de- 
pendence of molecular vibrational frequencies on bond 



9 



length was proposed empirically by Philip M. Morse in 
1929 [108] and later discussed by other authors [109]. 
We note, however, that there are significant differences 
between the systems considered here and those in the 
early work [108, 109] on the dependence of vibrational 
frequencies and the associated force constants on bond 
length: The vibrational modes considered in the previ- 
ous work[108, 109] were stretching modes, while mode 
IV considered here is torsional. Also the bond lengths 
considered previously[108, 109] were equilibrium bond 
lengths that differed because the bonds were in different 
molecules, whereas in the present work we are consider- 
ing non-equilibrium interatomic distances that differ due 
to differences in the geometrical constraints applied to 
the system. 

IV. CONCLUSIONS 

We have investigated theoretically the conforma- 
tions, the elastic low bias conductances and the inelas- 
tic tunneling spectra of molecular junctions containing 
propanedithiol and propanedithiolate molecules bridging 
gold electrodes. We found that, for junctions in which 
the molecule bonds to a gold electrode via a sulfur atom 
that retains its hydrogen atom, a relaxed structure of 
the junction could be obtained provided that the sul- 
fur atom bonds to a single gold atom of the electrode. 
This differs from propanedithiolate molecules bonding to 
gold via a sulfur atoms that have lost their thiol hydro- 
gen atoms in which case relaxed geometries with a sulfur 
atom bonded to two or three gold atoms have also been 
obtained [75, 78]. We have demonstrated that inelastic 
tunneling spectroscopy should be able to distinguish be- 
tween junctions in which neither sulfur atom retains its 
thiol hydrogen atom, one of the two sulfur atoms retains 
its thiol hydrogen atom, or both sulfur atoms retain their 
thiol hydrogen atoms. 

Comparison between our results for junctions relaxed 
with no geometrical constraints and the experimental 
inelastic tunneling data of Hihath et at. [26] indicates 
that most of their experimental measurements were car- 
ried out on junctions in which neither sulfur atom of 
the PDT molecule carrying the current retained its thiol 
hydrogen atom. We found no evidence of features in 
their experimental inelastic tunneling histogram [26] that 
might signal conduction via a propanedithiol molecule 
retaining both of its thiol hydrogen atoms. However, 



some weak features in the experimental inelastic tun- 
neling histogram [26] are consistent with the possibility 
of a small fraction of their junctions being bridged by 
molecules with one sulfur atom retaining its hydrogen 
atom, although molecules with gauche geometries and 
no thiol hydrogen atoms may also be responsible for 
these spectral features. Further experimental studies 
with greater sensitivity to weak features of the inelas- 
tic tunneling spectra are needed to distinguish between 
these possibilities. 

We also studied the evolution of the inelastic tunneling 
spectrum as a gold STM tip approaches an intact S-H 
group at one end of a propanedithiol(ate) molecule that 
is bonded to a second gold electrode via a sulfur atom 
that has lost its thiol hydrogen atom. We found that the 
frequency of a vibrational mode associated with torsional 
motion of the thiol hydrogen atom increases by approx- 
imately a factor of two as the STM tip approaches the 
molecule. Since we predict this mode to have a moder- 
ately strong intensity in the inelastic tunneling spectrum 
this frequency shift can be used to monitor the separa- 
tion between the gold STM tip and the thiol group exper- 
imentally as the tip approaches the molecule. When the 
thiol hydrogen atom detaches from the molecule we pre- 
dict this vibrational mode to disappear from the inelastic 
tunneling spectrum and the energy of the strongest IETS 
feature associated with vibrations of the sulfur to simul- 
taneously increase by ~4-5meV while its intensity also 
increases. 

Experiments observing this predicted behavior would 
be of interest since they would shed much needed light on 
how and when the thiol hydrogen atom detaches from the 
sulfur atom during bond formation between molecules 
and gold electrodes in single-molecular junctions. They 
may also help clarify the reasons why apparently similar 
experiments on gold-thiol molecular junctions often yield 
measured conductances differing by orders of magnitude 
[HO]- 
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